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ABSTRACT 

We calculate precise stellar radii and surface gravities from the asteroseismic analy- 
sis of over 500 solar-type pulsating stars observed by the Kepler space telescope. These 
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physical stellar properties are compared with those given in the Kepler Input Catalog 
(KIC), determined from ground-based multi-colour photometry. For the stars in our 
sample, we find general agreement but we detect an average overestimation bias of 
0.23 dex in the KIC determination of log(g) for stars with log(p)Kic > 4.0 dex, and a 
resultant underestimation bias of up to 50 % in the KIC radii estimates for stars with 
^kic < 2 Rq . Part of the difference may arise from selection bias in the asteroseismic 
sample; nevertheless, this result implies there may be fewer stars characterised in the 
KIC with R ~ 1 Rq than is suggested by the physical properties in the KIC. Further- 
more, if the radius estimates are taken from the KIC for these affected stars and then 
used to calculate the size of transiting planets, a similar underestimation bias may be 
applied to the planetary radii. 

Subject headings: stars: fundamental parameters — stars: oscillations — stars: interiors 



Introduction 



The primary task of the NASA Kepler Mission is to detect the perio dic dips in brightness 
caused by extrasolar planets transiting their host stars (jBorucki et al.ll2010l ). To accomplish this, 
the Kepler space telescope simultaneously monitors the brightness of up to 150,000 stars. While 
most of these stars are observed with a mean 'long' cadence of 29.4 minutes, a subset of up to 512 
stars at any one time may be observed with a mean 'short' cadence of 58.85 seconds. These high- 
precis ion, short-cadence observations make th e data suitable for asteroseismic studies of solar-type 
stars (jGilliland et al.l 120101 : iGarcia et al.ll201ll ). which can precisely determine the stellar structure 
from their resonant modes of oscillation. 

During the complete Kepler Mission, approximately 160,000 stars will be surveyed in a field 
covering over 100 square degrees. The Stellar Classification Project (SCP) was initiated in the 
pre-launch phase of the Kepler Mission to carry out ground-based photometric classification of th e 



targets in the Kepler field and to construct the Kepler Input Catalog (KIC) (IBrown et al.l 12011 



The primary task of the SCP was to accurately distinguish giants from dwarfs in order to select 
solar-type stars which may have detectable transits of planets orbiting in the habitable zone. The 
KIC also provided estimates of T e ff , log(g), [Fe/H] and E{B — V) from the analysis of multi-colour 
photometry, and derived properties (radius, mass and luminosity) from a simplified mass-radius 
relationship. Since most of the stars in the Kepler field-of-view had not been studied prior to the 
SCP, the classifications given in the KIC are often the only physical data available. 

In order to characterise the structure of a supe r-Earth, the planetary radius should be known to 
a precision of ~ 5% and the mass to within 10% ([Valencia et al.1 120071 ). To confirm the detection 
of an Earth-sized exoplanet based on its size requires the planetary radius to be similarly well- 
defined. As the transit signature only gives the planetary radius relative to that of the host star, 
it is essential to obtain an estimate of the stellar radius to a similar precision in order to fully 



-3 - 



characterise an exoplanetary system. 



The survey phase of the Kepler asteroseismology program took place during the first 10 months 
of science operations. Over this time more than 2600 stars were observed at short cadence, most for a 
duration of approximately one month. A small number were observed for the complete survey phase, 
giving higher frequency resolution in their acoustic power spectra. The asteroseismic properties 
of some of these stars have been analyse d in detail, identifying individual p-mode frequencies that 



2011: Mathur et al. 2011 



have been used to model t heir structures (jChaplin et al.ll2010i : iMetcalfe et alJl201Cl : ICampante et al 



The average asteroseismic parameters have been measured in a large 
ensemble of stars observed during the su rvey and have been u sed to determine physical properties 
and to test synthetic stellar populations ( Chaplin et al.ll201ll ). 



In this paper we calculate precise radii and surface gravities from the average asteroseismic 
parameters of 514 solar-type stars observed during the survey. These stars cover a range in apparent 
magnitude from 7.0 to 12.0 and a range in T e fj from 4400 K to 6900 K. We use these asteroseismic 
structural parameters to verify the accuracy and precision of the physical properties given in the 
KIC. 



2. Kepler Input Catalog 



There is insufficient telemetry bandwidth to monitor all of the stars in the Kepler field-of-view 
so targets are selected in order to maximise the probability of detecting a planetary transit around 
a solar-type star. Of particular interest is the detection of an Earth-sized planet in the habitable 
zone of a Sun-like star. Therefore, the primary task of the KIC is to correctly distinguish giants 
from dwarfs. 

The aim s and methods used for stellar classification in the KIC are described in detail by 
rown et al.l f)201ll ) . Here we give a brief summary. The physical data were derived from ground- 
based photometric observations in five bands (g, r, i, z, D51). The stellar values of T e g, log(g), 
[Fe/H] and E(B — V) were determined using Bayesian posterior probability maximisation to match 
the observed colours to stellar atmosphere models. The values of log(<7) derived from this method 
were then used directly to identify giants, defined as having log(g) < 3.6 dex. 

The KIC also gives estimates of radii, which were derive d from log (q) assum ing a mass-radius 



relationship that depends only on T e g . As stated explicitly by iBrown et al.1 (|201ll ) , this is true only 



in a statistical sense and takes no account of composition or evolutionary state. The accuracy of 
the primary stellar properties in the KIC is stated to be 0.4 dex for log (5) and ±200K for T e s, 
with the ability to distinguish main-se quence stars from giants estimated to be reliable to 98 % 
confidence for stars with T cff < 5400 K (|Koch et alJboid : brown et allboilh . 
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Asteroseismic data 



The frequency of maximum oscillation power, z/ max , and the average large frequency separa- 
tion, Az^, c an be obtained from the one-month Kepler ligh t curves using a number of different 



techniques (IMosser &; Appourchauxl 1 2009 ; 



2010; Mathur et al 



2010; 



Verner et al 



Huber et al 



2011 



2009 



Campante et al.l |2010| ; iHekker et al 



These characteristic asteroseismic parameters are 
sensitive to the globally-averaged stellar structure. To obtain estimates of lo g(q) and R from 



i/max and Au, w e used the Yale-Birmingham (YB) pipeline grid-based technique (jBasu et al 



Gai et al 



2011 



2010 



We tested the results obtaine d from the YB pipeline against those of similar 



methods (IStello et al.ll2009al : iQuirion et al.ll2010l ) and found them to be consistent 



The grid-based method determines the characteristics of stars by finding the maximum of the 
likelihood function of the set of input parameters {v max , Av, T e g, [Fe/H]} calculated with respect to 
a grid of stellar evolution models. The uncertainties on the derived log(g) and R were determined 
from the observational uncertainties on the input parameters and the priors that had been used to 
construct the stellar models. 

It should be noted that the seismic parameters calculated for the stellar models rely on the 
approximate scaling relations: 



M/M & (T eS /T eS ^ 3 - 5 



L/L e 



_l/ max,0) 







^ (M/M )°- 5 (T cff /r cff , ) 3 

A ^ A ^ 



(1) 

(2) 



Equation ([2]) assumes that the average large frequency separation scales as the square root 
of the mean stellar density. This relation has previously been tested using a wide range of stella r 



models (|Ulrichlll986l : IStello et al.ll2009bl : iBasu et al.ll2010l : IWhite et al.ll201lh . IWhite et al.l (1201 ll ) 



tested the validity of Equation (|2j) using models with a broad range of masses and evolutionary 
states and fo und a deviation fr om the scaling relation that is predominantly a function of T e g. 
According to IWhite et al.l (|201ll ) , any systematic error introduced as a result of the deviation from 
Equation ([2]) is largest for low-mass stars with T e fj 5600 K, where the value of Av calculated for 
the stellar models may be underestimated by up to 2.3 %, and smallest for stars with T e g- ~ 4900 K 
and T e Q ~ 6400 K, where any deviation from Equation ([2]) is negligible. For main-se quence stars 
with masses above ~ 1.2 Mq, which includes the majority of the stars in our sample. IWhite et al 



(2011 



state that it is best not to include their correction to the scaling law. As the stars in our 
sample cover a range of metallicities, masses and effective temperatures where any correction to 
the scaling relation is small, there should be no average systematic bias introduced as a result of 
assuming Equation ([2|). 

Equation ([1]) relies on the assumption that the frequency of maximum power should scale wit h 



the acoustic cut-off frequency (jBrown et al.lll99ll ; lKjeldsen &: Beddinelll995l ; lBeikacem et al.ll2011 



Unlike Equation ([2]), this has not been tested thoroughly with stellar models because there is less 
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confidence in the theoretic al predictions of the e xcitation and damping rates. The f ma x scaling has 
been tested empirically by IChaplin et al.l (|201ll ) who concluded that it should be reliable to a few 
percent and that the net bias on stellar properties as a result of assuming both scaling laws is no 
larger than the quoted uncertainties obtained from the grid method. The asteroseismic scaling laws 
ar e discussed in more detail and tested using Kepler observations of main-sequence and giant stars 



in lHuber et all ( 2011 



The KIC determinations of log (g), T e g and [Fe/H] were performed simultaneously, therefore 
to obtain a self-consistent comparison we used only the values of T e g and [Fe/H] given in the KIC 
as input to the seismic analysis. In the case of [Fe/H], the values given in the KIC are known 
to be imprecise. We therefore adopted a very large uncertainty of 0.5 dex in [Fe/H] so that the 
seismically determined values of log (g) and R were not significantly sensitive to metallicity. The 
large uncertainty in [Fe/H] is reflected in the derived uncertainties in log(p) and R. 

It has been suggested that the values of T e g given in the KIC are systematically too low for 
solar-type stars. Pinsonneault et al. (in preparation) have used an alternati ve analysis of the raw 
photometry in the KIC, based on the method described in I An et al.l (|2009l ). and found T e fj that 
are on average ~ 5% higher than those given in the KIC. A spectroscopic analysis of a sample of 
~ 100 stars has been performed by Bruntt et al. (in preparation) who found T e g values that are 
on average ~ 3% higher than those in the KIC. However, as it was our intention to test the KIC 
values of log(<7) and the radii derived from these, both of which depend on the KIC derived T e g, 
we used the KIC estimates of T e g to calculate the seismic values for log(g) and R. 

Using the KIC T e g provided a uniform set of temperatures for the stars under study and ensured 
that any systematic bias was not introduced as a result of using temperatures from another catalog. 
It should be noted that due to the weak dependence of log(g) on T e g, an increase of 5% in T e g 
corresponds to an increase of only 0.01 dex in the seismically determined log(g), and a decrease of 
only 1.2% in the seismically determined R. 

To test the robustness of the seismically determined values of log(g) and R, we applied the grid- 
b ased method to the sets of results for f ma x an d Av determined by five of the methods described 
in lVerner et al.l (|201ll ) that provided the greatest coverage over all stars with detected oscillations. 



The YB p ipeline was run w i th three different underlying mod el grids to determine the stellar 

giving up to 15 values of log(g) 



2008 



Gai et al 



2011 



properties (jDotter et al.ll200SI ; iMarigo et al 
and R for each star. 

The mean uncertainty in log(g) calculated using all methods and model grids was 0.025 dex. 
The mean standard deviation of the results from all methods and model grids was 0.007 dex. In the 
case of R, the mean relative uncertainty was 7.0%, while the mean standard deviation was 2.7%. 
The seismically determined stellar properties were therefore found to be insensitive to the choice 
of model grid or method used to determine ^ max and Av at the level of precision available. Since 
the scatter in the results was significantly smaller than the stated uncertainties, we compared the 
mean values of log(<7) and R to those in the KIC. 
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4. Verification of KIC stellar properties 

The distributions of the differences between seismic and KIC values of log(g) and R are shown 
in Fig. [TJ As the seismic radii span a range from 0.7 to 6.4 R@, the relative differences rather than 
the absolute differences are used. This necessarily causes a skewed distribution as the lower limit 
is fixed at —100%. The standard deviations of these distributions are 0.31 dex (log(g)) and 37% 
(R). The median values of the distributions indicate that there are significant biases in log(g)Kic 
(+0.17 dex) and R K ic (-18%). 

Since the radii given in the KIC were derived directly from log(g), an overestimation bias 
in log(g) causes an underestimation bias in R, since M = gR 2 (in solar units). Through this 
relationship, for a star with log(<7) = 4.0 dex, an overestimation of 0.17 dex in log(<7) results in an 
underestimation of 18 % in R. 



4.1. log(gr) 



A comparison of seismic and KIC values of log(g) is shown in Fig. [2J Assuming an uncertainty 
of 0.4 dex for log(g)Kic> we find H stars (2% of our sample) that would be identified as dwarfs 
using log(g)Kic hut are classified as giants using asteroseismology (these stars are plotted as crosses 
in Fig.©. 

We find that, on average, the KIC log(<?) follows the seismic log(g) when log(g)Kic < 4.0 dex. 
There is no bias apparent in the log(g) differ ences for these stars and the la agreement is at the 
0.27 dex level. This agrees with the results in lHekker et al.1 (|201ll ) who showed that the values of 
log(g) in the KIC are consistent with those determined from the seismology of large cohort of red 
giants observed by Kepler, with the majority having 2.3 < log(g) < 2.7 dex. The agreement becomes 
poorer for stars with log(g)Kic > 4.0 dex, with the values in the KIC on average being systematicall y 
higher than the seismically determined log(g). This behaviour was noted bv lBrown et al.l (|201ll ). 
who state that the KIC classifications tend to give log (g) too large for subgiants. 

There are two observational biases that may contribute to this discrepancy. Firstly, Malmquist 
bias increases the fraction of intrinsically brighter objects in the field, increasing the proportion of 
subgiants in the sample relative to their true number density. Secondly, the amplitude of oscilla- 
tions, and therefore their detectability, scales approximately with L/M, and therefore oscillations 
are preferentially detected for stars with lower log(g). These two observational effects combine to 
give a relatively high fraction of subgiants in the seismic sample. This is evident in the seismic 
log(g) distribution, which peaks around a value of 4.0 dex. The known shortfalls of the KIC lead 
to the overestimation of log(g) for these subgiants. The average overestimation bias for stars with 
log(g)Kic > 4.0 dex is 0.23 dex with a la scatter about this bias of 0.30 dex. 



Our findings confirm the warning given by iBrown et al.l (|201ll ) that users should be wary of 
log(g)Kic values for stars with photometric colours (g — r) < 0.65. Most of the stars in the seismic 



-7- 



sample fall into this range. The lower panel of Fig. [2] shows the difference between KIC and seismic 
log(g) as a function of (g — r). In the region where the KIC has shortfalls in the determination 
of log [g), the results from seismology can be used to complement the KIC results and perform an 
empirical correction. We note that the absolute photometry in the KIC is very reliable, with an 
estimated error of just 1.5% on the (g — r) colour estimates. 



4.2. Radius 

The KIC provides estimates of stellar radii derived from log(<7) and a simplified mass-radius re 



lationship. Asteroseismolo gy can provide more precise estimates of the stellar radius (jChristensen-Dalsgaard et al 



2010 



Metcalfe et al.ll2010l ). however most stars for which transits are observed will not have de- 



tectable oscillations. 

The relationships between log(g), mass and radius are shown in Fig. [3] for the KIC and seismic 
values. It is clear that there is more variation in radius for a given log (g) in the seismic data 
than appears in the KIC, which is not surprising because the KIC relationship is assumed to be 
valid only in a statistical sense. Through the KIC mass-radius relationship, the overestimation bias 
identified in log(<?)Kic (Sect. 14. ip appears as an underestimation bias in -Rkic- 

Fig. |4] compares the KIC radii with those determined from seismology. An underestimation 
bias of up to 50 % is apparent for stars with -Rkic < 2 Rq . Using the KIC radii would imply 
that 20 % of the stars on which we have detected oscillations have radii within 0.2 Rq of the solar 
value. However, the results from seismology prove that fewer than 3 % of the stars with detected 
oscillations are within this range. It must be stressed that the same observational biases described 
in Sect. 14.11 apply equally to the results for stellar radii. The cohort of stars that has been used in 
this comparison is subject to selection effects as this only includes the relatively small fraction of 
stars observed by Kepler that have detected oscillations. This is biased towards stars with higher 
L/M and may not be representative of all stars in the KIC. 

For stars with -Rkic > 2 Rq , the standard deviation of the distribution of relative differences 
between KIC and seismic radii is 35 % and no average bias is apparent. However, for stars with 
-Rkic > 2.5 Rq, a dichotomy is present whereby cool stars (T e g < 5000 K) tend to fall above the 
agreement line and hot stars (T e g > 6500 K) fall below it. These stars correspond to the extremes 
of the {g — r) plot shown in Fig. [2j For hot stars with (g — r) < 0.26, the bias reverses sign and 
log(<?)Kic is underestimated by up to 0.25 dex, leading to an overestimation of -Rkic by up to 50 %. 
Any empirical correction to the data in the KIC should therefore be done as a function of (g — r) 
and not as a function of log(g)Kic or -Rkic- 
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5. Discussion 

Using precise values of log(g) determined from the asteroseismic analysis of a large ensemble 
of solar- like oscillators, we have found that in its primary role the KIC effectively distinguishes 
dwarfs from giants. However, at the level of precision available to the asteroseismic parameters, we 
find an overestimation bias for stars with log(g)Kic > 4.0 dex, although this bias may be overstated 
relative to the full catalog due to selection effects in the asteroseismic sample. 

For stars with log(<?)Kic < 4.0 dex, on average we find an unbiased agreement between the 
se ismic and KIC value s, with a la deviation of 0.27 dex. This precision is higher than was estimated 
in lBrown et al.l (|201ll ) (0.4 dex). However, due to the presence of a significant number of subgiants 
in our sample, and the known shortfalls of the KIC that result in an overestimation of log (5) for these 
stars, a representative uncertainty of 0.4 dex for log(g)Kic is nevertheless probably appropriate. 

The radii given in the KIC are affected by the overestimated values of log (g) and as a result 
may be underestimated by up to 50% for stars with -Rkic < 2Rq. For stars with Rkic larger 
than this, the agreement is better with a relative la deviation of 35%. This is important to take 
into account if the KIC radii are used to estimate the number of ~ 1 Rq stars in our sample. The 
results from asteroseismology show that a much smaller number of stars with detected oscillations 
have radii close to 1 Rq when compared with the radii given in the KIC. This bias may lead to 
an underestimation of the radii of exoplanets transitin g the affected stars if the KIC radius is used 



to calculate the planetary radius (jBorucki et al.ll201ll ). however the bias is strongest for subgiants 



and may be considered an upper limit for less evolved stars. 

Although the KIC was not designed primarily to give values of log(g) and R to high precision, 
we find that the catalog may be used to give an indication of these stellar properties once the 
limitations outlined above are taken into account. Much better estimates of log(g) and R are 
available from asteroseismology for stars with detected solar- like oscillations. 



Acknowledgements 

GAV, WJC, YE and IWR acknowledge the support of the UK Science and Technology Facilities 
Council. SH acknowledges support from the Netherlands Organisation for Scientific Research. 
Funding for the Kepler Mission is provided by NASA's Science Mission Directorate. The authors 
wish to thank the entire Kepler team, without whom these results would not be possible. We also 
thank all funding councils and agencies that have supported the activities of KASC Working Group 
1. 



- 9 - 



REFERENCES 

An, D., Pinsonneault, M. H., Masseron, T., Delahaye, F., Johnson, J. A., Terndrup, D. M., Beers, 
T. C, Ivans, I. I., & Ivezic, Z., 2009, ApJ, 700, 523 

Basu, S., Chaplin, W. J., & Elsworth, Y., 2010, ApJ, 710, 1596 

Belkacem, K., Goupil, M. J., Dupret, M. A., Samadi, R., Baudin, F., Noels, A., & Mosser, B., 
2011, A&A, 530, 142 

Borucki, W. J., et al., 2010, Science, 327, 977 

Borucki, W. J., et al., 2011, ApJ, submitted f larXiv:1102.054"Tj) 

Brown, T. M., Gilliland, R. L., Noyes, R. W., k Ramsey, L. W., 1991, ApJ, 368, 599 

Brown, T. M., Latham, D. W., Everett, M. E., & Esquerdo, G. A., 2011, AJ, submitted 
(jarXiv:1102.0342|) 

Campante, T. L., Karoff, C., Chaplin, W. J., Elsworth, Y. P., Handberg, R., & Hekker, S., 2010, 
MNRAS, 408, 542 

Campante, T. L., Handberg, R., Mathur, S., et al, 2011, A&A, submitted 
Chaplin, W. J., et al., 2010, ApJ, 713, 169 
Chaplin, W. J., et al., 2011, Science, 332, 213 

Christensen-Dalsgaard, J., Kjeldsen, H., Brown, T. M., Gilliland, R. L., Arentoft, T., Frandsen, S., 
Quirion, P.-O., Borucki, W. J., Koch, D., & Jenkins, J. M., 2010, ApJ, 713, 164 

Dotter, A., et al., 2008, ApJS, 178, 89 

Gai, N., Basu, S., Chaplin, W. J., & Elsworth, Y., 2011, ApJ, 730, 63 
Garcia, R. A., et al., 2011, MNRAS, 414, L6 
Gilliland, R. L., et al., 2010, ApJ, 713, 160 
Hekker, S., et al., 2010, MNRAS, 402, 2049 

Hekker, S., et al., 2011, MNRAS, 414, doi:10.1111/j. 1365-2966.2011. 18574.x 

Huber, D., Stello, D., Bedding, T. R., Chaplin, W. J., Arentoft, T., Quirion, P.-O., & Kjeldsen, 
H., 2009, CoAst, 160, 74 

Huber, D., et al., 2011, ApJ, submitted 

Kjeldsen, H., & Bedding, T. R., 1995, A&A, 293, 87 



- 10 - 



Koch, D. G., et al., 2010, ApJ, 713, L79 

Marigo, P., et al, 2008, A&A, 482, 883 

Mathur, S., et al, 2010, A&A, 511, A46 

Mathur, S., et al., 2011, ApJ, 733, 95 

Metcalfe, T. S., et al., 2010, ApJ, 723, 1583 

Mosser, B., & Appourchaux, T., 2009, A&A, 508, 877 

Quirion, P.-O., Christensen-Dalsgaard, J., & Arentoft, T., 2010, ApJ, 725, 2176 
Stello, D., et al., 2009a, ApJ, 700, 1589 

Stello, D., Chaplin, W. J., Basu, S., Elsworth, Y., & Bedding, T. R., 2009b, MNRAS, 400, 80 
Ulrich, R. K., 1986, ApJ, 306, 37 

Valencia, D., Sasselov, D. D., & O'Connell, R. J., 2007, ApJ, 665, 1413 
Verner, G. A., et al., 2011, MNRAS, in the press (larXiv:1105.057ip 

White, T. R., Bedding, T. R., Stello, D., Christensen-Dalsgaard, J., Huber, D., & Kjeldsen, H., 
2011, ApJ, submitted 



This preprint was prepared with the AAS IATgX macros v5.2. 



-11- 



IB 




-1 -0.5 0.5 1 -100 -50 50 100 

log(g) KIC - log(g) seis (dex) (R KIC - R seis ) / R seis (%) 

Fig. 1. — Distributions of the difference between KIC and seismic values of log (5) (left) and 
radius (right). Error bars show the standard deviation in each range when determined using all 
possible combinations of model grid and input data. Dashed lines show the median values of each 
distribution. 
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Fig. 2. — Comparison of KIC and seismic determination of log (5). Left: Dashed line shows line 
of equality, dotted lines show the cut-off point for giants, crosses show stars identified as giants by 
seismology but not by the KIC, error bars show 20-point averages (excluding misidentified giants). 
Right: log(g) difference as a function of the colour index (g — r) taken from the KIC photometry. 
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Fig. 3. — Radius-log (g) relationship using data from the KIC (left) and seismology (right). Dashed 
lines show lines of constant mass for 1 M & (bottom) and 2 M (top). Dotted lines show the cut-off 
point separating giants from dwarfs. 




Fig. 4. — Comparison of KIC and seismic determination of stellar radius. Dashed line shows line of 
equality, crosses show stars identified as giants by seismology but not by the KIC, error bars show 
20-point averages. 



